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Motivated by recent scanning tunneling microscope (STM) experiments on cobalt clusters ad- 
sorbed on single wall metallic nanotubes [Odom et al, Science 290, 1549 (2000)], we study theo- 
retically the size dependence of STM spectra and spin-flip scattering of electrons from finite size 
ferromagnetic clusters adsorbed on metallic surfaces. We study two models of nanometer size ferro- 
magnets: (i) An itinerant model with delocalized s, p and d electrons and (ii) a local moment model 
with both localized d- level spins and delocalized cluster electrons. The effective exchange coupling 
between the spin of the cluster and the conduction electrons of the metallic substrate depends on 

■ the specific details of the single particle density of states on the cluster. The calculated Kondo 
coupling is inversely proportional to the total spin of the ferromagnetic cluster in both models and 
thus the Kondo temperature is rapidly suppressed as the size of the cluster increases. Mesoscopic 
fluctuations in the charging energies and magnetization of nanoclusters can lead to large fluctuations 
in the Kondo temperatures and a very asymmetric voltage dependence of the STM spectra. We 
compare our results to the experiments. 

?-j ■ PACS numbers: 73.22.-f,72.10.Fk,73.63.Fg 

CO 

■ I. INTRODUCTION 

S: 

Kondo physics [fij has seen a revival in recent years due to an exciting series of new experiments on mesoscopic [|| |^] 
and nanoscale||, ||, ||, |7j systems that has enabled a more thorough and controlled study of the basic problem of 
a local moment interacting with a sea of conduction electrons. It is now possible, for example in quantum dots, 
to tune the Kondo temperature, Tk, with externally applied gate voltages |2[ ^| or, with the scanning tunneling 
microscope (STM), to study the Kondo effect at single magnetic impurities on metallic surfaces.^, || The superb 
\ spatial resolution of the STM permits unprecedented direct local spectroscopic detail of Kondo impurities. 

Recent low temperature(~ 4 K) STM experiments have probed the spectroscopy of isolated nanometer and sub- 
nanometer ferromagnetic (Co) clusters on metallic single walled carbon nanotubes. [|| The tunneling spectrum of these 
small ferromagnetic clusters exhibited several interesting features-most notably Kondo effect in sub-nanometer diam- 
eter clusters and well resolved discrete level spacing in nanometer size clusters. Both sub-nanometer and nanometer 
size clusters exhibited a Coulomb charging gap in the spectrum near zero bias. The Kondo effect in the sub-nanometer 
size clusters appeared in the spectrum as a Fano-like sharp peak [Hjwith a half-width ~ 15 meV {Tk ~ 80 K) at 4 
K. This peak was not present at 100 K, presumably due to thermal destruction of the Kondo resonance. 

The goal of this paper is to present a theoretical framework to study the spin-flip scattering of conduction electrons 
from finite size ferromagnetic clusters adsorbed on metallic surfaces. In particular, we are interested in understanding 
some of the size- dependent features observed in the tunneling spectrum. We focus on how the Kondo Effect, the level 
spacing and the level widths depend on the size of the cluster. 

Cluster ferromagnetism may be described by two alternate models: (i) In the experiments of Odom et al.ffy 
magnetism of the cluster is most likely itinerant in nature as bulk Co is an itinerant band ferromagnet. The basic 
physics of itinerant clusters is captured in the model of Ref s. ft2| (see Sec. II). (ii) In some other cases, however, 
including the case of many semiconducting ferromagnets Jl3[ and rare earth materials, ferromagnetism is better 
described in terms of local moments that couple ferromagnetically or antiferromagnetically to the cluster conduction 
electrons (holes). This motivated us to introduce another exactly solvable model, where d-electrons produce highly 
localized magnetic moments. While the excitation spectrum of the two models turns out to be very similar, the 
way they couple to the metallic substrate is rather different. The most appropriate model to use in a given physical 
situation depends on the material that composes the nanoparticle and whether it is more appropriate to think of the 
magnetism of the particle as due to localized moments or itinerant electrons that, although free to move about, are 
nevertheless polarized strongly enough to give a net spin to the cluster. 

We make several assumptions during our computations: We assume the clusters to have between, say, eight and 
forty atoms and that they are large enough that bulk properties, such as the magnetization per atom, are not changed 
substantially due to their finite size. Furthermore, we assume that the single particle states on the island are of 
extended character, and therefore the single particle level spacing decreases as ~ 1/Na with an increasing number of 
cluster atoms, Na- We also assume that the ferromagnetic spin splitting A s of the extended states on the island is 
much larger than the single particle level spacing and is approximately independent of the size of the cluster. Thus 
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the total spin of the cluster is rather large and is roughly proportional to ~ Na ^> 1. 

Depending on the strength of the electron tunneling between the metal and the cluster we may distinguish two 
regimes: weak and strong. For weak tunneling we derive an effective Kondo Hamiltonian in second order perturbation 
theory in the tunneling. In this regime we find that at very low energy scales the cluster behaves most similarly to a 
quantum dot with a single unpaired electron. However, ferromagnetic interactions induce a large spin on the cluster 
and at the same time reduce the value of the effective exchange coupling between the spin of the cluster and the 
electrons in the metal. 

The sign of the exchange coupling depends on the specific details of the band structure on the nanoparticle. For a 
Co cluster, in particular, the itinerant model gives a tendency for antiferromagnetic coupling, and possibly produces 
a Kondo effect. However, the Kondo temperature Tk, decreases very fast with increasing cluster sizes, and therefore 
the Kondo effect can be observed only for very small clusters. 

The effective couplings we obtain, J eff , turn out to be in general M x M matrices describing spin- flip scattering 
among the M conduction band orbital modes that couple to the cluster (M=4 for metallic nanotubes[|i"4|). Therefore 
we expect that a series of Kondo effects takes place at Kondo temperatures corresponding to the antiferromagnetic 
eigenvalues of J eS . The total spin of the impurity, St gets successively screened at each Kondo temperature as the 
temperature is lowered: St St — 1/2 — > St — 1 .... For a nanotube this compensation can never be complete ]l5[ for 
any spin greater than 2 since there are only four channels available for screening. |fl4f The degeneracy of the remaining 
unscreened spin will ultimately be lifted by magnetic anisotropy induced by spin-orbit coupling. [ ]16| 

To make closer contact with the experiments of Odom et we compute the STM tunneling spectrum using 

the itinerant model. The simultaneous observation of the Coulomb gap and the discrete levels, whose widths almost 
equal the level spacing, suggest that the clusters investigated are in the intermediate tunneling regime. Nevertheless, 
carrying out perturbative calculations in the weak tunneling limit, we are able to reproduce the essential features of 
the STM spectrum above the the Coulomb charging energy. We find that small fluctuations in the charging energies 
can give rise to very asymmetrical STM spectra. 

From the STM spectra of the clusters at higher voltages, we can crudely estimate the Kondo temperature, Tk, of a 
small Co cluster of the size that exhibited a Kondo effect in Ref. fj| However, our estimates tend to give a Tk that is too 
small. We find that typical mesoscopic fluctuations in the charging energies of a ferromagnetic cluster can considerably 
increase the Kondo temperature, however, they do not account for the difference between our theoretical estimates 
and the experimentally observed Tk- [fj Nevertheless, the fact that some small Co clusters exhibit Kondo effect while 
others do not, indicates that these mesoscopic fluctuations may indeed play an important role in determining Tk- 
Since Tk depends exponentially on J cff , relatively small changes in the latter can induce large variations of Tk- 



This paper is organized as follows. In Sec. II we introduce the itinerant model of nanoscale ferromagnets. In Sec. Ill 
we study the spin-flip scattering of conduction electrons from an itinerant cluster using second order perturbation 
theory in the cluster-substrate tunneling, and we derive a general expression for the matrix of Kondo couplings, J . 
In Sec. IV we make contact to the experiments of Odom et aZ. by calculating in a simple approximation the STM 
tunneling spectrum of a ferromagnetic Co cluster on a metallic substrate in the limit of weak tunneling. In Sec. M 
we introduce a new local moment model for ferromagnetism and show how to compute the effective interaction J 
in this case. In Sec. we discuss in detail the connection between our work and recent experiments. Finally, in 
Sec. VII we present our conclusions. 



II. ITINERANT FERROMAGNETIC CLUSTER MODEL 

An itinerant mean field model of a magnetic nanoparticle was first proposed by Canali and MacDonaldjllJ and 
later extended by Kleff et a/.p^| to explain the dense spectrum in tunneling measurements of nanometer size Co 
particles. The basic Hamiltonian is 

^cluster = Yl e A° C 3° ~ ' ^ + ~ ^ ^ 

3,° 

where cJ CT (cj a ) creates (destroys) an electron with spin a at the j th energy level of the cluster with kinetic energy Cj 
(See Fig. |l|). The first term in Eq. (|l|) represents the kinetic energy of the hybridized s,p and d-band electrons. The 
parameter U > is the effective ferromagnetic exchange interaction on the cluster, Na is the number of atoms that 
constitute the cluster, and 

^ = 2 ^ c ]^ a 'C ja > (2) 

j a, a' 
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is the total spin of the cluster. The second term in Eq. (|l]) corresponds to ferromagnetic exchange on the cluster: 
It gives rise to a spontaneous polarization of the cluster by making spin alignment of different levels energetically 
favorable. Finally the last term describes the charging of the cluster with Ec the Coulomb charging energy, N — 
a ""i CT t ne total number of electrons on the cluster, and n g a dimcnsionless gate voltage that determines the number 
of electrons in the ground state. |Tg|| 

Throughout this paper we are neglecting fluctuations in the level spacing. In a real nanoparticle, the level spacing 
between the j th and j + 1 th energy level can be written as 6j = (Sj) + rjj, where (Sj) is the level spacing corresponding 
to the bulk energy-dependent density of states, (Sj) ~ l/g(ej), and r\j ((rjj)=0) is a fluctuating part due to shape 
irregularities in the nanoparticle. In this paper we take rjj = 0, corresponding to an infinitely strong level repulsion. 

The advantage of this model is that S z , S 2 and rij are conserved quantities and therefore the Hamiltonian in 
Eq. (Q) can be diagonalized exactly. In the more general situation a magnetic anisotropy term should be added to 
the Hamiltonian. Nevertheless, for the experimentally investigated clusters this anisotropy is estimated to be much 
less then the width of the Kondo resonance observed, and therefore it can be neglected. |l9[] 

The ground state can be explicitly constructed: jij) 

\S ,S* = Sa)% = I[ct T n c ]j vac > ( 3 ) 
3=1 3=1 

with 5*0 and N the spin and particle number in the ground state. The remaining states within the ground state 
multiplet may be explicitly constructed using the lowering operator: 



' 5 °' S -»g^ (2W„-'i)! ' S '»' S '''"'' S °- S °>g- <"» 

In itinerant ferromagnets there is an approximate rigid band splitting between spin up and spin down electron 
density of states. p(J Throughout this paper, we will assume that our clusters are large enough that the band 
splitting energy, A s , is well approximated by the bulk value of the specific material we are considering. In our model 
A s can be defined as the energy difference between the highest occupied spin up level, ca, and the highest occupied 
spin down level, ej: 

A s = e A - ei , 

and from band structure calculations it is known to be typically of the order of a few electron volts p0[ (see Fig. [l]). 
It can be determined by demanding that the ground state of the ferromagnetic cluster be stable to fluctuations of 
energy level occupations with constant particle number, and is related to the interaction parameter U as[[H[ 

U = ^A s + d Q , (5) 

Do 

where do is a small quantity that scales as 1/Na, while the first term is roughly independent of the cluster size, since 
So ~ Na- Using Eq. (||), one can estimate the energy cost of a particle-hole excitation (see Fig. ^B,C). This turns 
out to be of the order of the level spacing (oc \/Na) (and not the exchange splitting A s ): 

8E(So ±l,N) = mm{E cxcitcd {S ±1,N)- E G } ~ 5 A , &i ■ (6) 

where Eg denotes the ground state energy, and 5a (Si) is the level spacing near e A (ej). If U is large enough to fully 
polarize the cluster, the energy scales outlined above may change. |21| 

The minimum cost of adding a particle or a hole to the cluster can be defined as, 

SE ± . a = min{E(N±l,So + a/2)- E G }, (7) 

where the minimum is over all possible excited states. To determine SE±^ we consider processes like those in Fig. ||. 
These energies can be estimated as 

/ 1 e_ q a \ 

(8) 



(9) 
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where the corrections scale as 0(1/N\), and -B^ denote the charging energies in the limit of vanishing ferromagnetic 
coupling and an infinitely dense single particle spectrum: 

E^=E c {\±{N-n g )). (12) 

The parameter e = (e A + e/)/2 in Eq. (^) can be absorbed into the definition of n gi and we will set it to zero in what 
follows. 

The last terms in Eqs. (|^-pT|), proportional to A a /So, are specific to the mean field model discussed here and have 
little physical meaning. Their effect can be fully taken into account by renormalizing the 'chemical potential' e and 
the mesoscopic parameter do, and they therefore do not modify our results. 

Using a stability analysis similar to the previous one it is a trivial matter to show that typically 

8E±, a ~ E c /2 . (13) 

For later purposes it is also useful to have excitation energies of particles and holes added to any level defined, 
6E± a (N ± 1,Sq + cr/2). For example, increasing the total spin by adding an electron to the j th unoccupied level 
rather than the one immediately above e A gives, 

5E> +A =e 3 -(e A + S A ) + SE +A . (14) 

Removing an electron from level j below ej so that the total spin increases gives 

5E j _ A = ej - tj + 6E- T . (15) 

It is important to point out that in our model all effects of mesoscopic fluctuations have been put into the quantities 
Eci do and n g . All other quantities of the itinerant model are specified by the bulk band structure and the number 
of atoms, N A , in a given nanocluster. 



III. WEAK TUNNELING BETWEEN A FERROMAGNETIC CLUSTER AND A METAL 



A. General Expression 

A finite size ferromagnetic cluster can scatter electrons between different orbital channels in the substrate and flip 
their spin while doing so. If the spin-orbit coupling is sufficiently weak then SU(2) symmetry in the spin sector 
implies that the low-energy effective interaction between the cluster spin and the electrons takes the form (apart from 
a potential scattering term): 

•^Kondo = 2 ^2 Jfl "S ■ a^ ka a aa >a vk 'a' (16) 

k,k' 

where J 111 ' is a Hermitian matrix of spin-flip exchange couplings among the various orbital channels, indexed by /i and 
v. (In a one-dimensional system, such as a carbon nanotube, /i labels left-going and right-going modes of different 
symmetries. In an isotropic two or three dimensional host (i labels angular momenta about the cluster.) The total 
spin of the cluster is S and the a^ ka {a^ku) are creation (annihilation) operators for conduction electrons of the metal 
in orbital channel /i, wavenumber fc and spin a. A finite size cluster in an isotropic host will generally scatter electrons 
among all orbital channels, but in practice J MI/ can be truncated to include a finite number of channels, the largest 
angular momentum being I ~ where L is the size of the cluster and Xf is the Fermi wavelength of the conduction 
electrons. 
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In the limit of weak tunneling, J^ v can be calculated in second order perturbation theory in the tunneling. The 
Hamiltonian we consider is H = H + V, where H = H meta ,\ + Cluster and where H mcta i = J2^.,k,o- e ^ka\ krT a^ka is 
the Hamiltonian of the free conduction electrons of the metal and 

v = E (vt c l a ^ + h - c -) > ( 17 ) 

where 

V^ k = J2 V n ^(R n )^ARn) , (18) 
ra=l 

and N p is the number of points of electrical contact between the cluster and the metal. The tpj(R n ) are the wave 
functions of the j th levels on the cluster at the points of contact R n and the tpn,k(Rn) are the wave functions of 
the metal at R n . The tunneling amplitude V n can be estimated with a knowledge of the wave functions near the 
cluster-metal contact points and the work functions of the cluster and the metal. 
We can determine J^ v by equating 

(/fej-HE ^f'' . (19) 

n 

where \i) = \Sq, S z )q \v, ki, a) and |/) = \Sq,S z )q \h, kt, er') denote the initial and final states with energies Ei and 
Ef, respectively, and |n) stands for all possible intermediate states with energy E n , Ho\n) — E n \n). Here \v, k, a) 
denotes the Fermi sea with an additional electron in the vXh channel with momentum k and spin a. 

The easiest way to determine J^ v is to focus on purely spin-flip process to which potential scattering does not 
contribute, and choose S z = S a ,<r =i and S z ' = S - =|, giving (f\Hf ondo \i) = ±y/2S^.P lv . The RHS of 
Eq. (|l^) is evaluated by summing over both states with an extra particle and states with an extra hole in the 
intermediate state (See Fig. ||). After a rather straightforward computation we find that J 11 " is given by the sum of 
three contributions: 

The couplings J*i v and Jg v describe the contributions from tunneling processes involving doubly occupied and empty 
single particle levels: 



V 

■> A 



V V 3 u ' kf *V j ' ki I I (21) 

-1/2^ v ySE.j+er-ej SE- tl + A s + ej - e j J 1 ' 

— — V vt kf *v^ ki ( — 

+ 1/2^ " \8E +A -e A+1 



^ = > ^ ^ I - -e, - ^.,.A,' ,, : .,J > <*> 

and the incoming and outgoing electrons were put at the Fermi energy. Since A s 3> o\E±. CT both contributions are 
ferromagnetic, i.e., all eigenvalues of + Jj?" are negative. To see this let us first neglect the /c-dependence of the 
V£' k, s. Then the matrices J£ v (a = e, d) can be expressed as a sum: Jg" = Pj h ' ■ Each of these terms is obviously 
negative semidefinite, i.e. for any vector a M the product a^P^ a v is negative or zero. As a consequence, Jg u and 
are also negative semidefinite, i.e., in a diagonal basis all their eigenvalues are negative or zero. This simple proof can 
readily be extended for k-dependent V^f fc 's. 

The physical reason that the doubly occupied and empty single particle states give rise to a ferromagnetic contri- 
bution to 3^ v is that in the intermediate state in second order perturbation theory an electron (hole) that hops onto 
the island with spin parallel to that of the cluster has smaller energy than than an electron (hole) with opposite spin 
orientation, due to the ferromagnetic exchange on the cluster. Therefore an electron (hole) on the substrate with 
spin parallel to that of the cluster can lower its kinetic energy more efficiently through hopping to the empty (doubly 
occupied) states. This is reflected by the smaller energy denominator in Eqs. ( ^l|) and (|2^). 

Singly occupied levels, on the other hand, give an antiferromagnetic contribution to J^": 

Jr = ir E v ^* v i' ki (w — - — r- + 1*—; ) - ( 23 ) 



6 



with J£ v having only non-negative eigenvalues. p3| This is a result of the Pauli principle, because only electrons (holes) 
in states with spin antiparallel to that of the cluster may hop onto singly occupied levels in second order perturbation 
theory, and thereby reduce their kinetic energy. 

Eqs. (pi]), (^2|), and ( ^3|) constitute one of the central results of the paper. They describe how J^ v depends on the 
density of states of the singly occupied levels on the cluster, the excitation energies to the N ± 1 manifold of states 
and the tunneling amplitudes to the various levels of the cluster. 

In general, there is a competition among the three terms of Eq. (|2C|), and thus the sign of eigenvalues of depends 
on the specific structure of the single particle density of states on the cluster, and other mesoscopic parameters such 
as SE±^ a . In the absence of ferromagnetic interactions on the cluster, [7 = 0, however, the ferromagnetic contributions 
J^ v and Jj" identically vanish, and the effective interaction is always antiferromagnetic, provided there is an odd 
number of electrons on the cluster. 

The itinerant model has been used independently in Ref. [2^ to describe spin S = 1 metallic nanoclusters with 
similar results. In that case the ferromagnetic interaction is weak, U < S, the island is far from a ferromagnetic 
unstability, and the S = 1 ground state of the island results rather from two single particle states of the island having 
energies accidentally closer than U < 5. Then the ferromagnetic contribution of the doubly occupied and unoccupied 
levels is usually negligible, and the effective exchange interaction is dominated by the antiferromagnetic contributions 
of the singly occupied levels. B3| 

The off-diagonal elements of J**" are a sum of random numbers since the tp* (Rn)iifi,k(Rn) of Eq. (|l^ ) have random 
phase for different j and /z's. Therefore the size of the off-diagonal elements of will be down by a factor of 

i . / ( 5a -*) compared to the diagonal ones with [i — v, and the matrix J 1 "* will be dominated by the diagonal 



2S V A = 

elements if the number of scattering channels M is much less than 25*0. 



B. Weak Tunneling at a Single Point of Contact 



It is instructive and also useful to specialize the results of Sec. Ill A to the important case of a ferromagnetic cluster 
that makes contact with a metallic substrate in only a single point. Then we can choose a basis where J^ v consists of a 
single non-zero element as all but one conduction channel will decouple from the cluster (that is, a linear combination 
of the host orbital modes can be found such that only one orbital mode in the new basis has non-zero amplitude at the 
impurity location). With the origin taken as the single point of contact, Vj ,k of Eq. ( |l8| ) reduces to Vj = Vp* (O)V'fc(O). 
Eq. (p3|), e.g., then becomes 



1 



1 



So \SE +d - e I+ i + ej 5E- t i + 6A-ej 



(24) 



In the limit A s ,dE± -a 3> S Eqs. (|22|), (pi[), and ( p3[ ) can be approximately evaluated by assuming that the tunneling 
matrix amplitudes do not vary too much from level to level, and replacing the sums by integrals. For J s , e.g., we 
obtain: 



Js 



W\ 2 h 

So 



1 



1 



(25) 



with g(e) the single particle density of states on the cluster, and (|V/| 2 )j the average hybridization strength. The 
ferromagnetic couplings and J e can be expressed by similar integrals. For Sq 3> 1/2 these integrals can be 
combined in a straightforward calculation where the integrand SE ^ + ^ is separated into the two regions £ 3> 8E± j(7 
and £ -C 5E± a and then approximated in each region. We obtain the following estimate for the effective exchange 
coupling valid when 5E± tCr 3> 5: 



J' 



eff 



(\Vj 



So 



P 



p(e/)ln 



SE. 



+,1 



8E- 



(26) 



This simple expression is one of the central results of this work. Eq. ( p6|) describes the dependence of J cff on the 
particular details of the density of states of an itinerant ferromagnetic nanocluster. The first term involves a principal 
value integral (denoted by P) over energy and clearly shows that whenever ej < £ < £a the contribution to J is 
antiferromagnetic; otherwise it is ferromagnetic, i.e. electron and hole excitations in the singly occupied states give 
rise to an antiferromagnetic coupling while electron (hole) excitations of the empty (doubly occupied) states give 
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rise to a ferromagnetic contribution. The last two terms of Eq. fl26[ ) are mesoscopic fluctuations that depend on the 
specific charging energies of the nanocluster and the density of states at the top of the minority and majority bands. 
Other mesoscopic fluctuations come from variations in the tunneling matrix elements and in the level spacing, which 
are ignored in our model. These mesoscopic corrections become more pronounced for smaller cluster sizes and lead 
to strong fluctuations around J cff . 

For Co, the single particle density of states has a maximum within the spin-polarized part of the spectrum, ei < £ < 
€a, resulting in a J cff that tends to be antiferromagnetic. Numerical evaluation of Eq. ( |26| ) using the actual density 
of states for Cobalt shows that the ferromagnetic contributions to J cff from Jd and J e reduce the dominant 
antiferromagnetic contribution from J s by roughly 50-60 percent. Fluctuations of 8E± M in small clusters by a factor 
of 2 can lead to sample to sample fluctuations of ~10% in J eff for Co. This ultimately leads to large fluctuations in 
the Kondo temperature, Tk, from cluster to cluster, since J cff appears in the exponent of the expression for Tk- 

Eq. ( pif ) also describes how J cff scales with cluster size. From the normalization of the cluster wave functions we find 
(l^jl) 2 ~ ~n~- Assuming that the magnetization of the cluster per atom is determined by microscopic mechanisms 
and thus U ~ A s does not depend on the cluster size, both So oc Na and g oc Na- This results in a fast decrease 
of J eS with increasing cluster size (spin), J cff oc oc and therefore an exponentially suppressed Tk- Note that 
this suppression in Tk depends only on the level structure of the ferromagnetic nanocluster, not on any interference 
effects which may come from several points of contact, [pij 

It is useful at this point to mention what would happen if U <C §i,a and we had an odd number of electrons on 
the cluster. There is then only one singly occupied level on the cluster, the sums in Eqs. (21) and (p2|) vanish, and 

the sum in Eq. ( |23| ) reduces to just two terms. As a result, for the spin 1/2 cluster, J®? 2 oc ^ ^ gr^ oc N i- a with a > 

A 

since Ec oc l/iVg. Thus, for sufficiently large Na, J for th e ferromagnetic cluster (U ^> 8i : a) will scale to zero 
faster with Na than for the non-magnetic cluster. In Sec. IV B and in the conclusions, we will comment on what this 
means for the scaling of the Kondo temperature with the number of atoms in such a cluster. 



C. Weak Tunneling at Several Points of Contact 



For tunneling at several points of contact we again turn to Eqs. (20-23|). In this case, is necessarily a matrix 
reflecting the scattering of electrons in more than one orbital channel. As Eq. ( |l8| ) shows, there will be significant 
random fluctuations in the matrix elements J^ v from the ipj(R n ) as cluster size is changed. For N p points of contact 
there are a maximum of N p orbital channels that will be scattered. However, if N p is large, the number of orbital 
channels that are scattered may be much smaller. The largest angular momentum (orbital) channel scattered being 
I ~ j^, the matrix is approximately of size ~ I 2 x I 2 for a cluster in bulk or ~ I x I for a cluster in contact with 
a two-dimensional electron gas. For a nanotube at most 4 orbital channels can scatter low energy electrons. 

As we discussed earlier in this case, in principle at least, one may observe a series of Kondo effects. The cluster 
displays an under-screened Kondo effect, where the spin of the impurity is large and several channels with different 
coupling constants try to screen it. The channel with the largest coupling to the impurity screens half a spin first and 
ceases to interact with the impurity. Then the channel with the largest coupling from the remaining channels screens 
half a spin and so on. In practice, however, the various Kondo temperatures are exponentially separated, and the very 
small Kondo temperatures are likely to be much smaller then the spin-anisotropy energy of the cluster. Therefore in 
a realistic situation one can only observe one, or possibly two Kondo effects. p5[ 



IV. STM TUNNELING SPECTRA AND KONDO EFFECT 

A. STM Tunneling Spectra of Ferromagnetic Clusters 

In the simplest model of STM spectra the tip is treated as a point source of electrons and the resulting signal is 
proportional to the local density of states (LDOS) at the tunneling position r,]2^] 

^-{r,V) oc LDOS(r,u = eV) = -~-Im\G R (r,u = eV)]. (27) 

dV 7T 

Here G R (r, uj) is the retarded Green's function of the substrate below the tip, and V is the voltage difference between 
the STM tip and the substrate. Below we consider only tunneling directly into the cluster, ignoring any Fano-type 
interference effects. || [l(| 
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Let us first focus on the high-energy part of the spectrum at w > Tk, where the Kondo-type strong correlations 
are negligible. Using standard manipulations we can express the retarded Green's function as 



G G,E,E' \ \ 



1 



(w - (E - E G )) - M R (u) 



I (E'\tf(r)\G) 

EE' 



- ME) L + ( E -E a)) -M»J <E,| * W|G> I <28) 



EE' 



where now \E) denotes excited states of the system in the absence of tunneling, i.e. \E) is a direct product of the 
excited states of the isolated cluster and the excited states of the isolated substrate and E = 8e,e'E. The operator 
ip^(r) {0{r)) creates (annihilates) an electron at position r on the cluster. We average over the ground state multiplet 
\G) with degeneracy Nq = 2 So + 1, and M R is the generalized retarded self-energy matrix. 

Assuming that Tq = 2irgo Y^n < niin{i5£ , ±. (T , 8a, 8i}, we can do perturbation theory in the tunneling. In this 
case the ground state \G) can be approximated in leading order by the product of the independent ground states of 
the cluster and the metal, and the summation over the excited states \E) turns into a sum over the states \E± a ) 
defined in Sec. II (see Fig. ^). With these assumptions the retarded self-energy M R (ui) is approximately diagonal in 
E and E' 

M r (lj) ee , «real part - iSee'^Y. S ^ ~ E n)\( n \^\ E )\ 2 > ( 29 ) 



where summation is carried out over all possible intermediate states \n) and V is given by Eq. ([17]). The real part 
of the self-energy produces only a slight shift of the levels, and it can be neglected. The imaginary part, on the other 
hand, describes the broadening of the spectrum. 

Within these approximations the STM spectrum can be expressed as: 

where PL a is the total decay rate of the excited state with excitation energy 8E± a . 

Eq.([30|) assumes that no two states of the cluster in the same spin sector are closer to each other than r . Otherwise 
they may produce Fano-type resonances due to the interference between electrons (or holes) that tunnel through 
these states. For our mean field models with maximum level repulsion, however, this condition is always satisfied if 
To < Si,Sa- Also, if the ground-state energy difference for the cluster with N and iV± 1 electrons on it (or with spin 
So and So ± 1) is too small compared to To the initial state will not have a definite electron number (or spin) and our 
computation breaks down. 

For 8E± a > 8 a.i we find that the decay rate a is typically larger than the single particle tunneling rate ~ T J . 
This can be understood as follows: Assuming that the STM-cluster voltage is already large enough to overcome the 
Coulomb blockade, an energy ~ Ec/2 is transfcred to the cluster as an electron or hole hops onto the cluster from 
the STM tip. The extra charge carrier then tunnels into the leads. However, for 8E± a > 8a, i the island can be 

left behind in a state that contains electron- hole excitations. Assuming for the sake of simplicity that T — To is 
independent of the single particle state index j we obtain the result that the width of a given level is proportional to 
the number of ways it can decay: T J ± a ps N^^^To, where N^^ 7 " 1 is the number of island states with TV electrons 
whose energy E(N) — Eq < SE± a and which are accessible by the addition (or subtraction) of just one electron from 

the state corresponding to SE J ± a . There is no simple expression for the -^ccay^' 13 ' ^ey depend sensitively on the 
level j, Na, Eq and the fluctuating mesoscopic parameters do and n gi and we have to determine them for each set of 
parameters and each level separately. 

Within the mean field model we are using and to leading order in the tunneling V the number of decay channels 
of the state with energy E 3 ± a does not depend on j. This is an artifact of the mean field approximation: In higher 
orders of the tunneling (virtual tunneling processes) or in the presence of electron-electron interactions on the island 
(which may also change level occupations) excited states with larger energies generally have a larger decay amplitude. 
Examples of decay processes for a given excited state are shown in Fig. ||. 

According to Eq. ( |30| ) , the weakly coupled cluster spectrum is simply a series of Lorentzian peaks centered at the 
particle and hole addition energies determined by Eq. (|l|), Eq. (||) and the actual band structure of Co |2(J with a 
weight modulated by a random amplitude |yj(f)| 2 - In our calculation (See Fig. pj) of the STM tunneling conductance 
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we ignore (i) the amplitude modulation, (ii) fluctuations in the level spacing of the cluster beyond that given by band 
structure calculations ]2(| (as we discussed in Sec. II) and (iii) fluctuations in the hybridization rj w To. These effects 
could be included with some effort, but they would not be expected to change our conclusions. 

In principle, contributions from phonon scattering or electron-electron interactions also contribute to the relaxation 
rate T J ± a ; however, in the experiments of Ref. these can probably be neglected with respect to the single particle 
relaxation channel considered here. 

To compute the tunneling spectrum of an actual cluster with Na atoms we generated a discrete set of levels tj with 
a level spacing corresponding to the single particle density of states in Co.[^0j As discussed above, the non-uniform 
density of states is important to obtain a predominantly antiferromagnetic coupling, but is also important to obtain a 
stable ground state with a partially polarized band. The STM spectrum for low bias voltages is governed by the level 
structure near 6a and ei since these are the only states probed in an STM measurement at low voltage bias. The 
majority and minority level spacing at these energies can be estimated from band structure calculations as 5,4=5.55 
eV/N A , S r =lA3 eV/N A and S = 0.855A A .||} For Co, A s = e A - e 7 w 2eV. 

For the sake of simplicity we set e = (ej + e^)/2 = in our computations, effectively absorbing it into n g . This 
corresponds to a specific choice of contact potential between the lead and the cluster, and does not influence the 
overall features of the spectrum. 

Eq. (|3^) assumes the cluster is in its ground state before the electron from the STM tunnels into it. Thus, in the 
numerical calculation it is important to ensure that the parameters of the model give a ground state with definite So 
and A. Using Eq. (|l|) and Eq. (||) and keeping track of the l/N A corrections which are all put into do, one can derive 
the conditions for stability of the ground state, which we assume to have N s « 1.71Na singly occupied levels. By 
considering the fluctuations in the ground state spin shown in Fig. |l|B and C one can show that to guarantee stability, 
do must satisfy the conditions: 

N *x < Na (x M f<m 

For Co this simplifies to —1.67^- < do < 3.21 85 ^ A+1 ■ Similar conditions can easily be derived for the quantities 
Eq and Eq by considering fluctuations to the A ± 1 manifold of states like those in Fig. ^ and requiring SE±_ a > 0. 
For our model applied to Co we obtain Eq > — 1.0|p and Eq > 2.48-^-. 

A calculation of the tunneling spectrum of a ferromagnetic cluster with N A = 32 atoms and a corresponding spin 
So = 27 is shown in Fig. ||. For the computations we used the parameters Eq = .OleV, Eq = .08eV and To = .02eV 
in both figures. Then once do is given all the SE J ± a are determined via relations similar to Eq. (^|). The quantity 

A^dccty^ is computed numerically. In Fig. |5|A d = .07eV, and the ground state is close to being unstable against the 
lowest lying state within the (A, So + 1) subspace, while in Fig. [|B do = —0.05 eV, and the ground state is close to 
the lowest lying state within the (A, So — 1) subspace. 

In the spectrum in Fig. there is a series of sharp peaks at negative bias (electrons are removed from the cluster) 
while the spectrum at positive bias (electrons are added to the cluster) is rather smooth. These features correspond 
very well with asymmetrical features seen in some spectra of Ref. 0. On the other hand, the lower spectra of Fig. ^3 
has less contrast between positive and negative biases. The only parameter that has been changed between the two 
spectra is d of Eq. (H). Different values of d lead to different values of the number of decay channels Nfjjj^y 
for positive and negative bias, and hence lead to different widths of the peaks in dl/dV. This calculation shows how 
mesoscopic fluctuations may make a significant qualitative difference in the spectra of two clusters with the same Na, 
Ec and n g . For real clusters mesoscopic fluctuations in Eq and n g are also expected for fixed A^ and these will also 
lead to qualitative changes in the spectra. 

A careful comparison of our calculated spectra corresponding to the first few charging peaks with the experimental 
spectra shows that at voltage biases up to « ±Ec our calculations are in agreement with experiment, however at 
larger biases, eV> Ec (when two additional electrons (or holes) could be added to the cluster at the same time), 
the experimental spectra show additional structure, not present in our calculations, presumably due to many-body 
and/or non-equilibrium excitations left out of our simple model. 

As we will see in the next section, for the cluster size of Fig. || the estimated Kondo temperature is significantly 
below the experimental temperature, and no Kondo peak appears at zero bias. 



B. Estimating the Kondo Temperature 

It is obvious that Eq. (^0|) does not include Kondo correlations and therefore does not produce a Kondo peak in 
the spectrum. However, it is generally true that when the Kondo effect is present, it produces a peak of width ~ Tk 



10 



at the Fermi level with relative weight ~ due to the approximate unitary scattering, provided the temperature is 
less than T K ,T < Tx.0 || 



In order to estimate Tjc and its dependence on cluster size, we will use the results of Sec. Ill B for a cluster 
contacting the substrate in only a single point of contact. There are several reasons why we believe this may be a 
reasonable assumption to make for Co nanoclusters on metallic nanotubes.0 First of all, the ratio of broadening T 
of energy levels on the cluster compared to their separation ~ 6a j is experimentally found to be almost independent 
of the cluster size. Since 8a. i ~ 1/A^t and r ~ Np/Na (Np is the number of points of contact), this indicates that 
the number of effective tunneling points is approximately constant for the clusters investigated in the measurements. 
In fact, since the diameter of the clusters (.5-1 nm) of Odom et al. ^ are comparable to the diameter of a nanotube 
there is appreciable curvature at the cluster-nanotube interface, and it is quite possible that the cluster touches the 
surface of the nanotube only at a few points. Since the tunneling amplitude is exponentially sensitive to the tunneling 
distance, it is also possible that only one or two of these points dominates the conductance between the cluster and 
the nanotube, resulting effectively in a coupling in the form of a single point of contact. 

The standard expression for the Kondo temperature in the Kondo model is given by|l[] 



T K = DyJj cS Qo e- 1,r 00 (32) 

where D is the bandwidth of the conduction electrons of the metallic host and go is the density of states of the host 
at the impurity site. For a ferromagnetic cluster, however, this formula is incorrect, since the Kondo Hamiltonian 
provides an appropriate description of the cluster dynamics only below the characteristic energy of inelastic excitations, 
E- me \ = rmnjiJc', 8a, i}- As shown for example in the two-level system model, |^ above this energy scale fluctuations 
to various excited states destroy the coherent spin processes leading to the Kondo effect, and in the renormalization 
group approach in the regime above -Ei ne i the Kondo coupling remains unrenormalized. Therefore, for a single point 
of contact we estimate the Kondo temperature as 



T K ~ min{£ c , S A ,i}V J cS 9oe~ 1/r eo , (33) 

where we replaced the bandwidth D of the conduction electrons by the energy gap for inelastic processes on the 
cluster, D — > rtnn{Ec , 8 a,i} ■ Eq. (|3^) is only of logarithmic accuracy. In general, Tk contains an overall prefactor 
that incorporates corrections from higher excited states as well as possible mesoscopic fluctuation effects. This 
prefactor is usually of the order of unity, however in some cases it can be quite large and considerably increase the 
Kondo temperature. p8| 

In principle, the quantities Ec and 8a, i can be obtained directly from the STM spectra, and J go can also be 
related to the spectra via the width of the energy levels of the cluster. Recalling T Q — 2ir\V : >\ 2 Qo we can write 



J cS n 



2irS 



3C 



A 8 g(fl<% / \i I SE +d \ 8E +A 

p 1 « K-^-o - g{ei)ln [sE^r g{eA)ln se-; 



(34) 



where we have used Eq. (gfl) and taken \V j '\ 2 -> (| V 3 \ 2 ). Recall that g(e) cx N A . 

As is evident from Eq. (|3(j), the actual level width of cluster excited states observed in experiment is r j = iVdecayFg, 
where A^ocay is the number of energetically allowed decays for a particular excited state in the N ± 1 manifold as 

described in Sec. IV A. According to our model STM spectra calculations (Sec. |v]A) (which computes N^^l numer- 



decay 

ically) iVdccay is typically 1-3 for clusters with 7-30 atoms. As the most favorable case for obtaining an experimentally 
consistent Tk, we take iVdccay=l in our estimates. 

To obtain an estimate of the Kondo temperature for Co atoms adsorbed on nanotubes we used the experimental 
data of Ref. [?]. For a cluster experimentally estimated to have Na — 8 atoms, the value T ~ 0.24eV and Ec ~ .36eV 
can be directly determined from the experimental STM spectrum of cluster. Unfortunately, the level spacing 8i,a 
cannot be determined directly from the STM spectra for these small clusters, but can we estimate them by rescaling 
the level spacing measured at larger clusters, giving Si ~ 0.24 eV. 

Assuming that the spin splitting takes its bulk value, A s ~ 2eV we can then estimate Tk using Eqs. ( |34] ) and 
(|33|). We numerically evaluated both the integral in Eq. (|2^) as well as the discrete sums in Eqs. ( pl[j2^ ) using the 
actual Co density of states [p(| (assuming a single point of electrical contact). Neglecting the mesoscopic fluctuations 
in 8E± i(J , the integral and the discrete sum were found to be within 10% of each other for Na =8, 16 and 32. The 
Kondo temperature estimated this way for Na — 8, Tk ~ 0.16 K turned out to be about a factor 500 smaller than the 
experimentally observed Kondo temperature, T^ p ~ 80 K. (We obtain a value of J eS g = .12 for N A =8. This would 
need to be increased by a factor of ~ 2.5 to reach agreement with experiment.) Mesoscopic fluctuations of E± & may 
increase (or decrease) J eS by ~ 10%. Furthermore, our integrals (and sums) were evaluated with the assumption that 
Vj is independent of energy. This is not strictly true and will lead to additional fluctuations in J eS . These mesoscopic 
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fluctuation effects typically change the value of J cS altogether by 10— 15%. It seems unlikely such fluctuations could 
bring the theoretical estimate of Tr- into the experimentally observed range. Thus, there is a discrepancy between 



experiment and theory. There may be several explanations for this disagreement that we discuss in detail in Sec. VI 
of the paper. 

Eq. (^), on the other hand, is in qualitative agreement with the experiments in that it predicts an extremely rapid 
decrease of Tk with increasing cluster size. Since for a single point of contact Tq ~ 1/Na, Sq ~ Na, and g(£) ~ Na, 
the dimensionless exchange coupling scales as J cS go ~ 1/Na- We also verified numerically that the scaling J off oc 
of Eq. (j26) is maintained for the discreet sum in Eqs. ( B2H2!j ), 

Note that if one were to assume multiple contacts between the cluster and substrate, while keeping the decay rate 
fixed, then one would obtain smaller values for the typical tunneling matrix element V^ ,k and, hence, a lower Kondo 
temperature. 



As we found in Sec. IIIB , the dimensionless coupling constant for a non-magnetic cluster, J°J 2 So ~ l/N A - a , so that 
the Kondo temperature is suppressed less rapidly with increasing Na- In this sense, the large spin of a ferromagnetic 
cluster does not "help" the Kondo effect in anyway. 

V. LOCAL MOMENT CLUSTERS 
A. Local moment mean field model 

In many magnetic materials it is more appropriate to think of a localized d- or f-level (as in the Anderson Model) 
than to think of strongly hybridized s, p and d bands. These local moments may couple to each-other ferromagnetically 
and produce ferromagnctism. 

For small enough magnetic clusters with large enough Curie temperature, at low enough temperature the local 
moments form a large and rigid ferromagnetic spin, Sd, that couples to the extended states (a 'conduction band') of 
s- and p-character. The simplest Hamiltonian that one can conceive to describe this situation reads 

^ j Eq 

^cluster = X! e j C )<? C i° + j\T^ c ' ^ d + ~~T^ ~ n ^ 2 ' ( 35 ) 

with S c and Sd the total spin of the extended states and the local moments, respectively. A justification of Eq. (|35|) 
is given in Appendix [A]. Similar to the itinerant model, the first term describes the kinetic energy of the extended 
electron states, and the third term accounts for the finite charging energy of the cluster. The second term of Eq. (|35| ) 
describes the exchange interaction between the local moments and the extended states and tends to polarize the latter. 
We assume in what follows that the total spin of the conduction electrons is much smaller than that of the localized 
electrons. The exchange J is typically antiferromagnetic, J > 0, and the conduction electrons are polarized opposite 
to d-electrons.pp| 

The local moment model has St, S t (the total spin and its z-component), {rij}, Sd, and S c as conserved quantum 
numbers, and can thus be diagonalized exactly. The ground state is given by 

\St, S T )s c Sd = (S c ,Sd;Sc,Sd\ST,S T )\S c ,Sd;Sc,Sd) N , (36) 

where \S C , S d ; S*, S^ N = \S d , S^)\S C , S*) N and S T = S d - S c if J > 0. The state \S C ,S C ) N can be computed from 
Eq. (||) with S c replacing Sq everywhere, and \S C ,S^) N is computed from Eq. (||) with S c replacing So and S* replacing 
S z . 

Stability of the ground state implies the relation 

(Od + 1 j 

where A s is now the band splitting of the conduction band. 

Considering the type of particle-hole excitations shown in Fig. [j] and using Eq. (|37]) , one finds that the excitation 
spectrum is very similar to that of the itinerant model. In particular, we find that for charge fluctuations 5E± j(T ~ 
Ec/2, and spin fluctuations have a gap ~ 5a, Si. 
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B. Computation of the exchange coupling 



The coupling constants J^ v of the local moment model depend on the sign of the exchange coupling, J, of Eq. (|3 

1 

2 



Let us first focus on the case J > and Sd > S c . For the local moment model, (f\H^ ondo \i) = \^/2St J^ v ■ To 



evaluate the RHS of Eq. (|T^) we must again evaluate contributions from the conduction electron states with double, 
single and no occupation: 

For the sake of simplicity, let us consider the contribution of singly occupied levels, J£ v . Various matrix elements of 
the type ^ +1 (St + 1/2, St — l/2|cti |jSt, St) n arise in course of the evaluation of J^ u , and in contrast to the itinerant 
cluster model, the intermediate states of the local moment model (with J > and Sd > S c ) have an increase in total 
spin on the cluster in this case. To evaluate f +1 (S T + 1/2, S T - l/2|c] ; |S T , S T ) N we first expand \S T , S T )s c ,s d and 
\S T + 1/2, S T - 1/2)^+1 s using Eq. © to obtain 



f +1 (S T + 1/2, S T - 1/2|c}J5 t , S t ) n = 

So 

^ (S c , Sd; S*, St — S*\ St, St) 

S~=-Sc 



(S c - 1/2, S d ; S z c - 1/2, S T - S z c \ S T + 1/2, S T - 1/2) 

f +1 (S c -l/2,S|-l/2| C ]J5 C ,^) W . (39) 

What remains to be computed in Eq. ( |39| ) is the matrix element f +1 (S c — 1/2, S z — l/2|ct, |5 C , S Z ) N . To determine 
this, we use the states of Eq. (^) with S c replacing Sq. The overlap is computed by first directly evaluating ^ +1 (5 C — 
1/2, 5 C - l/2\c] l \S c ,S c ) N and then applying the Wigner-Eckhart Theorem for general S z . This yields f +1 (S c - 

1/2, S* - l/2|ctj_|5 c , S Z ) N = v/^|^. This can then be substituted into Eq. || which finally gives 

f +1 (S T + 1/2, S T - l^lcJjSr.Sr)* 
/ <? + S 1 - 2 

El c c j a q . nz c qz I q a \ 
y — ^ {b c ,^d,J c ^T - b c \bT,o T ) 

(S c - 1/2, S d ; S z c - 1/2, S T - S Z C \S T + 1/2, S T - 1/2) (40) 

The results of the evaluation of all the matrix elements on the right hand side of Eq. ( |l9| ) for the local moment model 
as well as the expression equivalent to Eq. ( p0| ) are somewhat lengthy so we relegated them to Appendix [B|. We 
evaluated them numerically for St = Sd — S c T(antiferromagnetic J and Sd > S c ) and found that the sing le p article 
contribution to the final exchange coupling differs in an overall sign from the itinerant model result of Sec. fl\ and it 
is ferromagnetic. On the other hand, if J is ferromagnetic, then St = Sd + S c and the matrix elements above agree 



in sign with those of Sec. Ill 



This can qualitatively be understood as follows: Suppose the total spin, St, of the cluster points upward. Then 
by assumption Sd > S c the spin of the d-levels, Sd, also points up. If the internal interaction, J, between Sd and 
S c is antiferromagnetic, delocalized electrons of the cluster with spin down will partially screen the local spin, Sd, 
so that the singly occupied states tend to have have spin down. A substrate conduction electron that hops on a 
singly occupied state must have, therefore, spin up that is parallel to the total spin, resulting in a ferromagnetic 
contribution to the effective interaction between the total cluster spin and the substrate. On the other hand, due to 
the antiferromagnetic interaction with the local spin, hopping to empty states with spin down have an energy smaller 
than those with spin up (parallel to the total spin), and give rise to an antiferromagnetic contribution to the effective 
interaction between the cluster spin and the substrate. The case J < can be understood along the same lines. 

Similar to the itinerant case, the signs of and J£ v are always opposite to that of J^ v regardless of the sign of 
J. Therefore, there is in general a competition between these terms, and the sign of the final coupling depends on 
specific band structure features. 

For an antiferromagnetic coupling, J > 0, in the limit where A s ,Ec 3> S, and S c ,Sd,ST 3> 1 we can obtain the 
following simple estimate for a single point of contact, analogous to Eq. (pi[): 



Jeff (l^f) 



St 



p L (e- " g{ei)ln [seTi + g{eA)ln set; 



(41) 
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The sign of the effective coupling depends on the sign of S c — Sd (here given for S c — Sd < 0). For completeness, we 
also give the expression for J < 0, which does not depend on the relative size of S c and Sd (note changes in mesoscopic 
fluctuations and overall sign): 



J Qo 



St 



P 



(£- ez )(e A -0 



- £>(ez)H 



SE. 



+.1 



5E- 



+ .T 



5E- 



(42) 



These expressions are particularly interesting: However large the constituent spins (S c ,Sd) of the cluster are, the 
cluster may still have a large effective coupling J cff if the coupling J between Sd and S c is antiferromagnetic (J > 0) 
and the total spin is sufficiently small. In most cases S c < Sd, therefore J eS is ferromagnetic and no Kondo effect 
develops. For S c > Sd, however, J eS changes sign (still assuming J > 0) and becomes antiferromagnetic. In this case 
a Kondo effect occurs with an effective coupling proportional to ~ I/St- 

It should be pointed out that the results described above are valid only in the weak tunneling limit. In the strong 
tunneling limit, the relative signs of J and J are switched. For example, when the "internal" interaction J > and 
the cluster is in the strong tunneling regime one cannot distinguish the cluster wave functions from those of the host. 
Thus, the J eff would have the same sign as J. Therefore we expect in this case J e g to change sign as one gradually 
increases the tunneling between the cluster and the substrate. We have not studied in detail how this transition would 
occur. 



VI. DISCUSSION 

This work grew out of an effort to better understand the experiments of Odom et in which Kondo effect 

was observed at low temperatures ~5 K for sub-nanometer Co particles adsorbed on metallic carbon nanotubes. 
Our model correctly predicts a Kondo temperature that decreases quickly with increasing cluster size, however our 
numerical estimates of Tk tend to be too small by a factor of ~500 for a cluster with Na=8. 

There may be several explanations for our low estimate of Tk- 

(a) Mesoscopic fluctuations in Eq. (^6|) can eventually increase the effective Kondo coupling and thus bring Tk close 
to its experimental value. However, since the sign of the mesoscopic fluctuations is random, this interpretation would 
appear to contradict the experiments in which a significant fraction of small Co clusters produced a Kondo effect. 

In addition to fluctuations in the various charging energies, the tunneling parameters Vj and the level positions tj 
also fluctuate from cluster to cluster. These additional fluctuations were neg lected in Eq. <M), since their contributions 
decrease with increasing cluster size. For small clusters, however, they may produce important additional fluctuations 
inJ off . 

(b) It appears furthermore that the experiments were performed close to the mixed valence regime as the width of 
the levels is comparable to the Coulomb charging gap. The effective Kondo Hamiltonian we derived in second order 
perturbation theory may not adequately predict Tk in that case. In general, approaching the mixed valence regime the 
Kondo temperature becomes larger than expected by the naive Kondo model calculation, the Coulomb gap shrinks, 
and the Kondo resonance gradually merges with the high-energy part of the spectrum. 

(c) It could be possible that some of the Co atoms in the cluster are not strongly attached to the others, and in the 
STM spectrum one observes the signal of these individual atoms. This explanation is, however, very unlikely in our 
opinion, because the Co atoms show a very strong tendency to cluster formation, and moreover the Kondo resonance 
is observed rather uniformly over the surface of clusters which are supported on nanotubes. 

(d) In our analysis, we assumed that the anisotropy energy is smaller than the Kondo temperature, and therefore 
neglected it. However, it is conceivable that very small clusters have a considerably larger anisotropy than our 
estimates based on experiments on large clusters. [0 According to the STM measurements, (tJ Co clusters in the 
nanotube experiments tend to have 'pancake-shape' and the relative position of the Co atoms is probably strongly 
modified with respect to the bulk due to the presence of the substrate. Although the value of spin-orbit interaction on 
Cobalt is not particularly large, it is still possible, that the highly anisotropic shape of the cluster and the deformed 
bonds generate an anisotropy, that is larger or comparable to the observed Kondo temperature, Tk ~ 70K. The 
effect of anisotropy on the behavior of the grain is rather complex, and we shall discuss it in a subsequent publication 
p4j . We would like to mention, however, two important results that may be relevant to the experiments. Large 
spin anisotropy is usually unfavorable to the Kondo effect: In most cases it leads to an Ising-like behavior with 
exponentially suppressed effective Kondo couplings, and gives rise to a dramatic decrease of Tk- However, for very 
small grains with a half-integer total spin and an almost perfect planar anisotropy, it can result in an effective strongly 
anisotropic Kondo coupling that is considerably larger than the couplings in Eq. (^4|) . We find that for the smallest 
grains in Ref. |?] strong planar anisotropy could give rise to a Tk in the experimental range. 

(e) Another possible source of error is our assumption that the calculated bulk density of states can be used for a 
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small cluster. If the peak in the density of singly occupied states is shifted significantly from the energy value shown 
in Ref. 20, the value of J might be increased. It is also possible that many-body corrections, omitted from our mean 
field model, could increase the value of J cS sufficiently to account for the discrepancy with the experimental Tk ■ 

It is interesting to compare the results for ferromagnetic Co clusters with nonmagnetic Ag clusters studied on single 
wall metallic nanotubes. [?J The Ag clusters showed no Coulomb gap or discrete level spacing in the STM spectrum. 
This suggests that the Ag clusters were not in the weak tunneling regime where valence fluctuations can be ignored, 
and where an effective Kondo Hamiltonian can be derived for a particle with an odd number of electrons. (Indeed, 
if the coupling to the substrate is sufficiently large, the mean number electrons on the cluster may be far from an 
integer, and the distinction between even and odd becomes meaningless.) Our analysis suggests that a Kondo effect 
can occur for a particle of a non-magnetic metal, with odd electron number, if the coupling to the nanotube is in an 
appropriate intermediate regime. 

It is also interesting to compare the results for Co particles on nanotubes with measurements of several Co particles 
on a highly-oriented pyrolytic graphite (HOPG) sheet reported in Ref. |7. The STM measurements did not show 
apparent single particle levels in the latter case. Assuming that the coupling to a nanotube and graphite were not 
too different, this could be explained by a higher density of states on the HOPG surface, |n) (Recall that for 1 nm Co 
clusters on nanotubes, the level broadening was roughly equal to the level spacing. Thus the levels may be broadened 
beyond resolution on the HOPG surface.) The STM measurements typically show a minimum in the dl/dV spectrum 
near zero bias, when tunneling into the Co cluster on HOPG, but the width of the feature is relatively large. When fit 
to a Fano formula for a Kondo resonance, the authors of Ref. obtained values of Tx of order 700 K even for clusters 
as large as 1 nm in diameter. Since it was not possible to raise the temperature enough to observe a temperature effect 
on the tunneling feature, however, supporting evidence for existence of a Kondo effect could not be obtained from 
this source. We note that STM measurements for tunneling directly into the HOPG substrate also show a mininum 
at zero bias. 

Differences in the coupling of Co clusters to a nanotube or graphite surface may also play a role in the observed 
spectral differences. Theoretical and experimental studies of STM images of graphite surfaces have indicated that 
there is an asymmetry in the local density of states at nearest neighbor atoms. |32f] This difference may also play a 
role in the interpretation of the spectra of Co on HOPG. Finally, it is possible that the matrix element for coupling 
between the cluster and the nanotube is reduced relative to the coupling to graphite due to the curvature of the 
nanotube. 



VII. CONCLUSIONS 



In this paper we have studied electron scattering from ferromagnetic clusters on a metallic substrate. We studied 
two cluster models. The first model describes itinerant ferromagnetism |ll[ and is probably appropriate for the 
description of experiments such as those of Odom et aZ.|?J on Co clusters. We also proposed another solvable cluster 
model, where spins on the d-levels are treated as localized entities. This latter model may be more appropriate for 
nanoscale rare earth ferromagnets or semiconducting ferromagnets such as GaMnAs, though in both cases spin-orbit 
interaction plays an important role and leads to strong spin-anisotropy effects. 

We derived a general expression for the Kondo couplings J^ v for both ferromagnetic cluster models. The sign of 
the obtained coupling depends in both models on the details of the band structure. For the itinerant model, virtual 
tunneling onto the singly occupied levels on the cluster induces an antiferromagnetic exchange interaction, while 
doubly occupied and unoccupied levels generate a ferromagnetic contribution to the exchange coupling. 

We have shown that for Co clusters the itinerant model leads to dominantly antiferromagnetic coupling between 
the cluster spin and the conduction electron spins. However, fluctuations to doubly occupied and empty states give a 
large ferromagnetic contribution to the exchange coupling that reduce it to roughly half its original value, and thus 
cannot be neglected. (As we discussed in Sec. [II, for the non- ferromagnetic spin 5=1 islands studied in Ref. 
these ferromagnetic contributions are small.) 

The exchange coupling J^ v involves various scattering channels. Therefore, in principle, the cluster could produce 
a series of Kondo effects where the spin of the cluster is gradually screened. p5[ 

It is important to emphasize that in the regime of weak electron tunneling between the metallic substrat e and the 
cluster, ferromagnetism has no special role in producing the Kondo effect as we have emphasized in Sec. IV B. In 
fact, our calculation shows that ferromagnetism tends to suppress the Kondo temperature with increasing Na more 
so than for the case of a non-ferromagnetic cluster. Besides the "strength" of the ferromagnetism, U, the Kondo 
scale is also affected by the density of states on the cluster and finite charging energy, as well as the cluster-metal 
conductance. 

The weak tunneling analysis we performed is only appropriate if the conductance between the cluster and the metal 
lead is smaller than the quantum conductance. Increasing the number of tunneling points leads to an increase in 
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the cluster-metal conductance. Once this conductance becomes larger than the quantum conductance, the effective 
charging energy is renormalizcd to a value close to zero, and a perturbative computation in V breaks down. In this 
regime extended states on the cluster are strongly hybridized with those in the metal, and can be viewed as part of 
the extended states in the metal. 

In the regime of strong electron tunneling between the substrate and cluster, it is not clear whether the itinerant 
model is able to produce a Kondo effect. On the other hand, our local moment cluster model gives a natural description 
of the strong tunneling regime. In the local moment model the localized d-electron spins can be viewed as a magnetic 
cluster embedded in the metallic host. This model has been analyzed in detail in Ref. ^4] In our local moment model, 
antiferromagnetic exchange (J > in Eq. ([35])) between the local moments and the conduction electrons produces 
a Kondo effect in the strong tunneling regime, though, the Kondo temperature decreases very fast with increasing 
cluster size. We have also argued that within the local moment model, the effective coupling between the electrons 
in the substrate and the cluster spin must change sign as one gradually increases the tunneling between the cluster 
and the substrate. 

Both the itinerant and local moment calculations show that the Kondo coupling is inversely proportional to the 
total spin St of a ferromagnetic cluster, which in turn, is proportional to the size of the cluster. |35|1 The Kondo effect 
is due to quantum fluctuations of the cluster spin, and these are suppressed as 1/ St for large ferromagnetic clusters. 
Thus Tk goes to zero exponentially with increasing cluster size. 

To make stronger contact with the experiments of Ref. |^ on Co clusters on a carbon nanotube, we also calculated 
the STM spectra of a ferromagnetic cluster as described in Sec. [i"v|A. We found that mesoscopic fluctuations in 
the charging energies may give rise to interesting qualitative changes and asymmetries in the STM spectrum. It is 
possible, for example, that the positively and negatively charged states of the cluster have very different decay rates 
and therefore the positive (negative) voltage side of the spectrum shows discrete levels while the negative (positive) 
voltage side displays a continuum spectrum. 

We can use the model parameters extracted from the high energy part of the STM spectra to make an estimate 
of Tk and predict how it scales with cluster size. Our results agree with the experiments in that they produce a 
rapid decrease of Tk thereby rendering the Kondo effect impossible to observe in larger clusters. However, the Kondo 
temperature we find is already too small by a factor of ~ 500 compared to the Tk observed for a small cluster of ~ 8 
atoms. In Sec. VI we have enumerated a number of effects which might raise Tk relative to the predictions of our 
simple model which may be a possible explanation for the discrepancy between theory and experiment. We remark 
that by neglecting the ferromagnetic contributions to the exchange coupling one would obtain a Tk that is larger by a 
factor of ~ 10 3 . Thus the ferromagnetic contributions to the effective coupling are essential and cannot be neglected. 

Many open questions remain regarding the physics of small ferromagnetic clusters. Among them are: (i) Accurate 
estimates of the net spin of 5-50 atom clusters supported on a substrate, (ii) Magnetic anistropy energies in clusters 
this size, (iii) The nature of non-equilibrium and other many body effects. We believe that the STM is a crucial tool 
for gathering cluster-specific data for ferromagnetic nanoparticles and will undoubtedly reveal even more intriguing 
physics of these tiny systems in the years to come. 
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APPENDIX A: DERIVATION OF THE LOCAL MOMENT MODEL 

Consider Na magnetic impurities embedded in close proximity in a metallic host. The Hamiltonian is H = 
-Hmetai + -f^int where -ff me tal = Ylu k a ^.fc'Vfco- describes the free conduction electrons and 

H in t =-J2 Jf ^ ?%{r) ■ S d (r>) + J ]T ■ Sd(r) (Al) 

r,r' f 

describes the direct interactions the impurity d-levels and the interactions between the conduction electrons and the 
magnetic impurities. The first term in Eq. ( |Al| ) describes the ferromagnetic interaction among the localized d-lcvels 
of the impurity atoms and the second term, Hk = J^2pS c (r) ■ Sd{r), describes Kondo scattering from these d-levels 
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by the conduction electrons. Jjr(f, r") > is the ferromagnetic exchange interaction between two localized d- levels 
and J > is the bare Kondo exchange coupling between the d-levels and the conduction electrons. The conduction 
electron spin operator at position r is 



9 H <P*j (r)c] a a aa 'c ra , 



(A2) 



where ipj(f) is the wave function of conduction electrons with level index j at r. The utility of Eq. (Al) is two- fold: 
(i) It gives an expression for the important limiting case of a "cluster" which consists of just one impurity. In this 
limit, = JS c (0) ■ Sd(0) which is just Hk for a single impurity, (ii) Eq. (Al) can describe the limit of a cluster so 



strongly coupled to the metallic host that the conduction band conduction electrons of the cluster and those of the 
host cannot be distinguished. Thus, the character of the <Pj{f) that appear in Eq ( |A2| ) vary depending on the physical 
situation. For a single impurity, ifj if) is the wave function of the host metal conduction electrons, for a cluster in the 
weak tunneling regime as described in Sec. Ill <Pj(r) describes the conduction band wave functions of the ferromagnet 
or ferromagnetic semiconductor and for a cluster in the strong tunneling regime ifj (r ) is a strong hybridization of the 
ferromagnetic conduction band electrons and the conduction electrons of the host metal. 

To derive Eq. (35), we neglect anisotropy in the cluster magnetization formed by the localized d-levels. We further 
neglect spin-wave excitations to states of different total d-electron spin of the cluster. In a Hciscnberg-type model one 
may estimate the cost for such an excitation as -E S pin-wavc ~ Jir(afe) where a is the lattice spacing of the cluster's 
atoms and k is the largest wave vector of the spin wave allowed by the physical size of the cluster. If we take Jf ~ 



0.1 eV, e.g., then for a small cluster of size L ~ 10 Awith lattice constant a ~ 2.5 A, k n 



so that the minimum 



spin-wave energy £™£_ wave ~ 60meV » 
can concentrate on the subspace where 



600K. Below that energy scale spin waves can be therefore neglected and we 



Y / S d (r) = N A S d = ST 



(A3) 



Within the subspace of maximum total d-level spin, the d-level spin is a rigid spin which can only change its projection 
on the z-axis. The collective effect of the impurity d-levels is to give the cluster a net spin. It is this spin that conduction 
electrons will scatter from-either "directly" in the limit of strong tunneling between the cluster and metallic host or 
"indirectly" as described in Sec. Ill in the limit of weak tunneling. (In the limit of weak tunneling only conduction 
band electrons may hop on and off the cluster.) Within the S d = S'™ ax subspace, 



(S d ,S'£\H K \S d) S* d ) = ~ J2 ^M^4a ff aa>C fa (S d ,S'£\§ d ^\S d ,S^ 



r,j,j' 



9 [r)c] a a aa ,Cj, a ^-{S dl S' d z \S d \S d , 



r,j,j' 



The sum over r can be estimated, 

f 

Neglecting the off-diagonal terms, 



j = j' (normalization) 
j j' (random numbers) 



J 



s c ■ s d 



The full Hamiltonian is 

-^cluster 



E 



N A 



S c ■ s d 



E c 



(£'• 



^2j F (r,r)S d (r) ■ S^) 



(A4) 
(A5) 

(A6) 

(A7) 
(A8) 

(A9) 



where we have put in the Coulomb charging energy by hand. Recall that now the c\„ refer to conduction band 
electrons, not s, p and d hybridized bands as in the itinerant model. The last term in |A9| is just an irrelevant shift in 
the total energy in the subspace S d = 5'™ ax so we drop it. Thus, we arrive at 



^cluster — X e j%<r + ~j\j~Sc ' S d H ^"(£ fyc 



(A10) 
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APPENDIX B: MATRIX ELEMENTS AND KONDO COUPLINGS FOR THE LOCAL MOMENT MODEL 



In this Appendix, we include some more lengthy expressions not incl uded in Sec. VB. To complete the evaluation 
the RHS of Eq. ( JTgj ) for the local moment model described in Sec. [VB we need the matrix elements: 



M 2 



N (S T , S t - % T |S T + 1/2, S T - l/2)f +1 

Sc I s — S z 
- E V °2S C (S c ,Sd; S*, St - S* ~ 1\St,St - 1} 

sj=— s c 

{S c - 1/2, S d ; S z c + 1/2, S T ~S Z C -\\S T + 1/2, S T - 1/2) 



(Bl) 



Ms 



f-\S T + 1/2, 5 T - l/2|c jT |5 T , St)^ 



So 

£ 



25, 



—(S c , Sd; S z , St — 5*1 St, St) 



(S c - 1/2, S d ; S z c - 1/2, S T ~ S Z C \S T + 1/2, S T - 1/2) 



(B2) 



Ma 



N 



N-l 



(S T , S T - M^St + 1/2, S T - 1/2) j 
/ , \/ C oc c (S c ,S d ;S z ,ST - S z - 1\S T ,S T - 1) 



s;=-s c 



2S C 



<S C - 1/2, S d ; S z + 1/2, S T - S* - 1\S T + 1/2, S T - 1/2). 



(B3) 



The matrix elements then directly yield the expression for the generalized Kondo couplings for the local moment 
model: 



V 3=1+1 \ 



M 4 M 3 



M 2 M X 



(B4) 



where Mi denotes the matrix element already given in Eq. (^0|). The products M4M3 and M2M1 can be evaluated 
directly: 



jr = - 



St 



1 A 



1 



1 



j=j+i 



(B5) 



giving a ferromagnetic contribution when J > 0. If J < and St = S^ + Sc (T - >!) one finds — s 1 +1 
agrees with Eq. (p3|) of the itinerant model with St — > So- 



A calculation similar to the one that led to Eq. (B4) and Eq. (B5) then yields: 



■J- which 



Jr = Y, v t ,kf * v » M] 

j>A 

and from the doubly occupied states 



2St + 1 



2S C 

S T + 1 V 2S C + 1 



SE +:i - e A +i + ej 8E+,t + A s - e A+1 +tj )' 



1 ( 2S C \ 
St + 1 \ 2S c + 1 J 



JUis _ \ ~> V J,kf* v j,ki 2S T + 1 



(B6) 



(B7) 



where we have agained assumed J > 0, giving > and j% v > since 8E± :(T ~ Ec/2. In the limit of a single 
point of contact, one recovers an expression similar to Eq. (26), Eq. (|4l|), except with an overall sign difference when 
J > (and the precise form of the mesoscopic fluctuations). In the case J < (take T~^l an d |— >•?) wc have 

2S7+T "> -237+T and "> which makes J e v < and < 0. 

The reduced matrix elements |33|] for the all the states of the local moment model are given below: 
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1. J>0 



?-\S T + l/2\\ Cj \\S T )Z = 



2S T 



V2S t + 1 



L (5 T + l/2|| C ,||5 T )f = y/2S T + 2 
^ T + l/2|| c j||S T }f = 



V2S T + 2 V 2S C + 1 
f +1 (S T -l/2|y||S T )f = v^TT 



2. J<0 



V2SVTT 

?-i(s T -i/2\\ Cj \\s T )% = Vm^J^l- 



f + i<S T -l/2||cj||S T >f = ^§+i 



f- 1 ^ - l/2|| Cj -||^>f = 
f +1 (5 T + l/2|| C ]||5 T )f = v^T+l 

l (^-l/2|| c t||5 T )^ 
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FIG. 1: Spin excitations of a ferromagnetic nanoparticle. For the precise definition of €aii, Sa/i, an d A s see text. Fig. A: The 
fully porarized ground state. Figs. B and C: Lowest lying particle-hole excitations having energy ~ Sa/i- 
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FIG. 2: Charging excitations of a ferromagnetic particle. Circles indicate particles added to (removed from) the ground state. 
Fig. A (B): Lowest energy state for adding a majority (minority) electron to the cluster. Fig. C (D): Lowest energy state for 
adding a majority (minority) hole to the cluster. 
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FIG. 3: Kondo scattering from a ferromagnetic cluster. Figs. A and C: Time ordered diagrams in second order perturbation 
generating the single particle contribution Eq. ( |23| ) to the effective exchange J M ". Continuous lines represent the incoming 
and outgoing electrons while dashed lines denote the intermediate excited state of the nanoparticle. In A an electron in state 
\v, k, l> hops on the cluster onto level j and another electron leaves the cluster from level j with outgoing state \fi, k' , f>. In C 
first an electron hops out from level j into state k' , j> and then the incoming electron hops onto level j. Level occupation 
changes are shown in D. 
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FIG. 4: The dashed lines represent possible possible decay channels for the states \E+^){N + 1 electrons) and |-E-,t) (N — 1 
electrons) to lower energy states with N electrons. In an actual decay, only one of processes indicated by the dashed lines 
would occur. 
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FIG. 5: Calculated STM spectra of nanometer-size ferromagnetic clusters with parameters Na = 32, Eq = .OleV, Eq = .08eV, 
To = .02eV, and fi = at 4 K (see main text). We took F = 27r|Vj| 2 go ~ To, independent of j, and \ip(f)\ 2 = const in Eq. (fjo|). 
In Fig. A we set do = -07eV while in Fig. B do — — -05eV. The difference between the two cases demonstrates the sensitivity 
to mesoscopic parameters. 



